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lave been recently found to exhibit many t( 
complex networks. It has been reported tl 
3wed and directed. In this paper, we report c 
like model when the topological features oi 
i of real biological networks. Specifically, ini 
^, we deal with a complex heterogeneous net 
gree distribution coupled to a continuous dy 
rhe dynamics of the model is very rich and 

are observed upon variation of the model's 
3S numerically and report on several quanti 
L and size distributions of clusters of stationa 
[Its are discussed in view of recent debate ab( 

nature and on the basis of several biologica 
IV the dvnamics studied. 



and the world- wide- web, the challenges i 
cation strategies in order to provide fas1 
I] , the implementation of better algorithn 
,tion dissemination [12,13], and the under 
ith the final goal of protecting the networ 
tial attacks and virus spreading [6,14,15,1 

!scribed by several characteristics. Among 
3d, one usually finds that real-world netw 
aeans that the average distance between 
i^ith the system size only logarithmically 
xity of these networks is characterized by 
connectivity k) of a given element (node 
world networks studied so far display a 
3wer-law P{k) ~ k~^ (termed scale-free r 



; all levels of organization are nowadays th 
d theoretical research. Recent analysis of ] 
las provided new useful insights into biolc 
janization [19]. It is also believed that a 
Drotein networks will help to elucidate tl 
)teins whose functions are unknown [20]. 
discovery of how biological entities inter 
cesses and tasks, as well as how their fu 



3 steady, periodic and chaotic states tha 
m's parameters. The results obtained all 
s about the robustness (hereafter intended 
d the phase space of chaotic dynamics) ■ 
iological networks. 

r is organized as follows. In section 2, w 
a, introduce the model and explain the ni 
;nt dynamical regimes are shown and anal 
3tion rounds off the paper by discussing ( 
i of the present study. 



icuss in what follows is built in two layers 
■ of the underlying network while the sec 
^namics of the network's components. Ai 
levant biological networks has been recent 
'his characteristic is shared by other netw 
are directed. Henceforth, we assume that 
rk corresponds to a biological entity and 
^tions. 



; and discuss the results for two plausit 
5sion and reaction kinetics. 

tivity of the nodes is described by the 
i}, where Qi {i = 1, . . . , N) accounts for th 
ide i in a network made up of N elem( 
described by the following set of first-or^ 



fF(G(t)), 



ne nonlinear term where the interactior 
•e taken into account. Equation (1) inclu 
boolean Networks [33,34] as well as conth 
ks [35], both widely used to model perio( 
logically relevant situations. Additionally, 
s-Menten description [3,32,36], 

action matrix introduced before. Additioi 
Cj is the connectivity of node i, and the fi 



)ut to be very rich and, depending on the 
iptotic dynamical regimes are observed, cl 
Qd chaotic attractors. All three states m? 
dization, each in different islands or clust' 
interconnected through nodes which have 
iptotically) their dynamics are effectively 

periodic states point to regions of the p; 
letworks might operate, the existence of c 
al, inconsistent with the reproducibility c 
organisms. Hence, we have characterized 
1 under variations of both h and p and 

probabilities of ending up in either chac 
he distribution of clusters or islands of nc 
Dver, due to recent interest in what is knc 
aalyzed the topological features of the clut 
ior. 

isented in the rest of the paper were devel 



f Qi are taken from a uniform distribution 

F the equations is performed using a 4*^ 
The total integration time is large com] 
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having chaotic, Pch or periodic Pper behavi 
meter p. The value of h has been set to 4 ai 
jdes. See the text for further details on the < 



tworks of sizes ranging from N = 100 to i 
id h, we have performed at least 1000 it 
I time step in the integration scheme wae 
Dn a further criterion in order to obtain tl: 
iodic states. 
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lies of the nodes is not uniform across the 
iitj in the initial conditions and that of 
3 nodes reach a stationary state, others 
;ic behavior. The following argument exj 
a node i is such that Wij = — 1 for all j, 
The same will happen for those nodes / 
• ?''s of the previous kind, etc. . . Now, tv, 
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IS that the uniform state of node inactivity b' 
a function of the parameter p. 10^ different r 
le of p. 

state yielding a positive value of A. On t 
loes not end up in a chaotic regime and ( 
ds are observed. Consequently, the prob 
stained, but periodic orbits are observed 
total number of realizations in which A < 
3rbit. 

probabilities as a function of the topolo^ 
h = 4 and A^ = 300. Two different threshc 
ily determined by Pper, Pi ~ 0.2 and p2 ~ 
londs to the case in which most of the i: 
dividual dynamics are described by frozei 
: 0). On the other hand, when the inters 
;ory {p > P2), the activity of the nodes c 
ilq. (1). In the intermediate region, pi < p 
xj, periodic and chaotic) are achieved. 
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of the quantity gi[gi = 0] (activity level whei 
node with the parameter h. This node belc 
ifurcations when increasing the value of h s1 
1 of the whole network at /i = 0. In the 
of nodes in a stationary state with non-zei 
^6 the nodes of this cluster end in a periodi( 
Lyapunov exponent of the network as h is i: 
I two bifurcations, (c) Periodic trajectory h 
)onding to the node of figure (a), the value 
3vel of the same node to the stable fixed po 
bifurcation) when the initial condition of th 
1 in (c). 

we h = h aX> which the largest eigenvalue 
;tors such that all their components have 
5ji + 5hWij becomes positive. Then h is 
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mch of stationary solutions exchanges stc 
3a). The computed largest Lyapunov exf 
Fig. 3b near h ~ 0.345: it approaches zero 
Aon parameter value, and then decreases 
ongs to the new stable stationary brand 
lay non-zero activity. As shown in Fig. 3a, 
increases with h. Eventually, this state be( 
typically through a Hopf bifurcation (ei 
rtate in which the activities oscillate (Fi^ 



both Pper and Pch in figure 1 are zero outs 
3 parameter p clearly indicates that non-st 
interplay between excitatory and inhibitc 
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.t a node with periodic dynamics converges 
(in arbitrary units). The results are average 
id at least 100 different initial conditions f 
= 0.7. h has been fixed to 4. 

confirmed by computing the frequency 
:s converge to a periodic cycle for differe: 
lure is as follows: First, we identify the 
punov exponent is zero. Then, we focus o: 
e identified, a vector Tl^ = {t\,f2, . . . ,f^} 
)eriodic dynamics gi. The t*'s stand for tl 
a) = 9i{ti) = ■■■= gi{ti) and dgi{t\)/dt -- 
[n this way, after verifying that t* — i!_i i 
)nding i-orbit is given by this constant. 

le probability that a periodic cycle has ; 
a in the figure, it is very likely that the ft 
island lies around ujper = 1. It is also c 
3tric, but biased towards larger frequenc 
m explanation to this behavior. It may pi 
stribution of the nodes and the specific v; 
lumber of input and output connections ; 
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a node ending up in a chaotic state as tl 
leriod doublings starting from a periodic c 
)f h (from (a) to (f)) are: 5.30, 5.50, 5.63, 
A^ork parameters are as in Fig. 6. 

ues in the phase diagram for the chaotic r 
!.. Clearly, as the degree of nonlinearity inc 
•e frequently, which translates in a largei 



;riodic cycle, which doubles its period su( 
phase. This corroborates that when h a: 
the behavior of the system is dominated 
or chaotic). Moreover, the fact that X) Pc 
ivhich are not in a chaotic state may be 
this regime of parameters, when a given 
is very likely that it has periodic clusters 

.escribed the activity patterns in terms ol 
3r, one of the most interesting aspects of r 
topology is highly heterogeneous. There c 
lany others. This should have some beai 
next section, we round off our numerical 
the dynamical properties unraveled with 
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;e the structural properties of distinct dynj 
^nitudes. The first and simplest structure 
istribution of nodes exhibiting either peri 
tograms of periodic (and chaotic) cluster 
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t a connected cluster of nodes displaying e 
given size (in number of nodes forming the ch 
to two radically different behaviors, they a 
e as in Fig. 5. 

ident dynamics appears as one of the mos 
. In Fig. 9, we show four different clusters 
dynamics. As one can see, no typical stri 
mparable sizes, except that all of them hi 
ster average connectivity. 

i like to add a few sentences about net 
)me of utmost interest in biological and i 
LS in modularity and community structui 
)onents or loops that appear more frequ( 
. identical degree distribution [3,42]. We 
le clusters of periodic or chaotic behavior 
e underlying network. In particular, chac 
' a universal topological pattern, contrai 
3r studies on network evolution and mot 
in ref. [3], and [43]). We believe that thi 
ing network, although heterogeneous anc 
lom scale-free network, for which the pro 
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laving periodic or chaotic behavior in a node 
hat dynamical regimes are well differentiated 
leters are as of Fig. 6. 

e have monitored the probabiUty that a : 
ither chaotic or periodic behavior. Thesi 
espectively, are defined as the ratio betwe 
k' greater than or equal to k that end i 
i the total number of nodes with conne 
ealizations. Figure 10 shows the results o 



terns have been studied in the last several 
2,3,6]. This approach has allowed the un 
topologies in many well-studied problems, 
ith other topologies and those for real grap 
g of epidemic diseases [15,16] or rumor 
networks to random failures and attack; 
r plays a fundamental role. 

ms, a very rich behavioral repertoire is w( 
:al systems range from circadian clocks to 
itration of certain chemicals, for instance, 
ilysis. On the other hand, the molecular b 
en discussed, though chaotic behavior is 
)gy as intrinsic noise makes it difficult tc 
current experimental techniques. The dy: 
could plausible describe at least to biolo 
on and reaction kinetics in metabolic net 

ne dynamics has to rely necessarily on n 
n to that of real genes. However, there £ 
the individual dynamics of gene expres 
Fman several decades ago, one of the mos 
Qdom Boolean Networks (RBNs) [31], w] 
le predictable properties and guided our 
x descriptions. Our model belongs to a g( 



;ic attractors is likely. By studying simpl 
here — the intrinsic complexity of the pn 
and detailed description of real gene dy 
if the parameter space (i.e. {p,h)) that b 

on this hypothesis, by exploring magniti 
, we would either guess dynamical inter. 

experimental validation of the structui 
tter seems more feasible due to latest di 
es, biocomputational tools, and data colle 

e, the continuous model employed here si 
those seen in RBNs. For instance, althoi 
3f h entering Eq. (2), it lies in the interv 
stence of a relevant average input conne 
logeneous networks in the context of RB] 
by Kauffman [31], who later suggested 1 
(up to (k) = 5) if the Boolean functions 
homogeneity [45]. Our results provide ai 
igh value of (k). Namely, that heterogei 
ighly nonlinear interactions allow for lar 
) what has been observed in RBNs, pi{h 
L size, which makes the previous analysis 



works, the system of differential equatio 



se experimental results are only now beco 
;e several features of interest such as that 
parameter range in which the existence 
at distinct behaviors correlate with conm 
jld be tested when more experimental ii 
study may help understand biological pr( 
eaction kinetics with the tools of network 
J^inally, we would like to point out that f 
ither random scale-free networks with low 
ous nets) [46], suggest that the main con 
tly depend on 7. 
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